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Abstract 



We study a Yukawa theory with spontaneous chiral symmetry breaking and with a large 
number N of fermions near the finite temperature phase transition. Critical properties in 
such a system can be described by the mean field theory very close to the transition point. 
We show that the width of the region where non-trivial critical behavior sets in is suppressed 
by a certain power of 1/N. Our Monte Carlo simulations confirm these analytical results. 
We discuss implications for the chiral phase transition in QCD. 



1 



1 Introduction 



The transition in QCD separating the high temperature quark-gluon plasma phase from 
the low temperature hadronic phase has been studied intensively in the last decade. Un- 
derstanding the properties of this transition is becoming increasingly important in view of 
recent experimental progress in the physics of ultrarelativistic heavy ion collisions. 

Since the u and d quark masses are small the dynamics of the finite temperature transition 
is affected by the phenomenon of chiral symmetry restoration which occurs in the limit when 
the quark masses are put to zero. In this limit QCD with two massless quarks has a global 
SU(2)ixSU(2)ij symmetry which is spontaneously broken to SU(2)y at low temperatures. 
It can be argued that if the restoration of this spontaneously broken symmetry proceeds 
through a second-order phase transition the critical properties of this transition are in the 
universality class of the classical 0(4) spin model in three dimensions [fll 0]. This means that 
the leading singular behavior of thermodynamic quantities can be predicted, i.e., it is given 
by universal 0(4) critical exponents. However, universality does not answer more detailed 
questions such as how this criticality is approached and what is the width of the region of 
parameters in which this singular behavior sets in. These questions require more detailed 
knowledge of the dynamics of the theory. 

In this paper we discuss a phenomenon which is related to the way the critical behavior 
sets in for a certain class of theories with second-order chiral symmetry restoration transition 
at finite temperature, T c . These are theories with a large number of fermion species N. On 
general grounds, one could expect that the universal critical, or scaling, behavior near T c sets 
in as soon as the correlation length, £, of the fluctuations of the chiral condensate exceeds 
1/T C . However, as we show in this paper, if the number of fermions involved in the chiral 
symmetry breaking is large, the critical behavior which sets in when £ exceeds 1/T C is given 
by the mean field theory, rather than by the arguments based on dimensional reduction and 
universality as in The critical behavior given by these latter arguments sets in much 
later, closer to T c , when the correlation length £ exceeds N x /T c . Below we shall determine 
the value of the positive exponent x. 

The phenomenon of the non-trivial critical region suppression has been observed in the 
Yukawa model || and in the Gross-NeveuQ model using large- N expansion and was 
confirmed by lattice Monte Carlo calculations in [|J . These large- N results predicting mean- 
field critical behavior were in an apparent contradiction with the general arguments of [|TJ. 
In this paper we show that both critical regimes are realized in the vicinity of T c , but in 
separate scaling windows, one following the other. 

2 Large N Yukawa theory near T c 

In this section we consider a general Yukawa theory in d dimensions, 2 < d < 4, with 
a large number N of fermion species and at finite temperature T. As argued in ||, in the 
continuum limit both Yukawa model and Gross- Neveu, or Nambu-Jona-Lasinio, models with 

1 In this paper we shall use the terms "Gross- Neveu model" and "Nambu-Jona-Lasinio model" inter- 
changeably, especially when it concerns a theory with a four-fermion interaction in 2+1 dimensions. 
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a four-fermion interaction define the same theory. In the absence of a bare fermion mass 
there is a (chiral) symmetry, which can be broken spontaneously at low temperature with a 
suitable choice of couplings. This symmetry is restored at some finite temperature T c . We 
are interested in the nature of this phase transition. 

In the absence of the fermions (or when the Yukawa coupling is zero) the nature of the 
transition is rather well known. It depends on the symmetry of the model which is restored 
at T c . In this paper, for simplicity, we consider a model with Z 2 symmetry. Similar results 
will also apply to theories with other symmetry groups (e.g., SU(2)xSU(2), as in QCD with 
2 massless quarks), as long as the temperature driven symmetry restoration transition is of 
the second order. 

One can argue that the critical behavior of a quantum theory of a scalar field near T c is the 
same as in a classical scalar theory with the same symmetry. The argument is based on two 
expected properties of the model: dimensional reduction and universality. In the Euclidean 
formulation the quantum scalar field is defined in a <i-dimensional box with the extent in the 
imaginary (Matsubara) time dimension equal to 1/T. When the diverging correlation length, 
£, becomes much larger than 1/T C long wavelength fluctuations of the field become effectively 
(d — l)-dimensional, i.e., on their scale the box looks like a (d — l)-dimensional "pancake". 
Since such fluctuations determine the critical behavior in the model one can expect that, by 
universality, the critical exponents are the same as in a d — 1-dimensional theory with the 
same symmetry. This d — 1-dimensional theory is obviously a classical field theory at finite 
temperature. One can also understand this realizing that the classical thermal fluctuations 
whose energy is 0{ksT) dominate over the quantum fluctuations with energy O(hou) for soft 
modes of the field. 

Another common way of describing this phenomenon in perturbation theory is to consider 
Fourier decomposition of the field into discrete Matsubara frequency components. Non-zero 
frequency acts as a mass term of order irT for the (d — l)-dimensional components of the 
field. Near T c this mass is much larger than the mass of a component with zero Matsubara 
frequency. The dimensional reduction is then equivalent to the decoupling of the modes with 
nonzero frequencies. 

We want to understand what happens in this theory near T c when one turns on the 
Yukawa coupling. Dimensional reduction and universality arguments suggest that the critical 
behavior of the theory should not change. This is based on the observation that there is no 
zero Matsubara frequency for the fermion fields due to antiperiodic boundary condition in 
the Euclidean time. Therefore all fermion modes should decouple at T c . In other words, the 
fermion fields do not have a classical limit and do not survive quantum-to-classical reduction 
atT c §. 

However, as was demonstrated in J|, the fermions do affect the behavior near T c in a 
certain way. We shall show that this happens when there are "too many" of them. The 
theory can be solved in the limit when the number of fermions, N, is large. In this limit 
the theory has mean-field critical behavior near T c 0, [|. This is different from the critical 
behavior of the corresponding classical scalar field theory. The mean-field behavior was also 
observed in numerical Monte Carlo calculations near T c [|J. 

Here we show that such mean-field behavior can be reconciled with the standard ar- 
guments of dimensional reduction and universality. The phenomenon which leads to an 
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apparent contradiction is the suppression of the width of the non-mean-field critical region 
by a power of 1/N. 

We consider the following model with one-component scalar field and Z 2 symmetry in d 
dimensional Euclidean space: 

11 N - 

£ = ^? + + A0 4 + E + 9<P) 1>f (!) 

We regularize the model by some momentum cutoff, A. The cutoff can be removed if 
d < 4 but we shall keep it finite to compare with a corresponding lattice theory. There 
are two other important scales in the theory: the temperature, T, and the physical mass, 
m, which we identify with the mass of thermal excitations of the scalar field. Near T c this 
mass, m, is significantly different from the zero temperature mass mo, since m vanishes at 
the critical temperature. It is the mass m, or the correlation length 1/m, which is important 
for the critical behavior. The mass m is a natural measure (more natural than, say, T — T c ) 
of the distance from the criticality.0 Therefore, near the finite temperature phase transition 
we have the following hierarchy of scales: A ^> T ^> in. 

Let us consider the renormalization group (RG) evolution of the couplings from the scale 
of A down to the scale of T and then from T down to m. We focus on the quartic self coupling 
of the scalar field, A. The evolution from the scale A to the scale T is governed by the RG 
equations of the ^-dimensional quantum Yukawa model. After that, at the scale of T, we 
pass through a crossover region due to the fact that the fermions and nonzero Matsubara 
frequencies of the scalar fields do not contribute to the evolution below T (the decoupling). 
The evolution below T is governed by the RG equations of the scalar 4 theory in d — 1 
dimensions. 

If the window of scales between A and T is wide enough (as it is in the continuum limit 
A — > oo) the value of the renormalized coupling A at the scale T, A(T), is close to the infrared 
fixed point of the <i-dimensional Yukawa theory. In the large- N limit one can calculate this 
value 0j: 

C4 - d) T 4 ~ d 

A(T) ~ ^ for 2<rf<4. (2) 

The case d — 4 is special. The infrared fixed point is trivial and is approached logarithmically 
as A/T — > oo: 

This value provides the starting point, A^-i = TA(T), for the evolution of this coupling 
below T in the 4 theory in d — 1 dimensions. For large N this coupling is small. As we shall 
see shortly, this is the reason why the critical region where one observes non-trivial critical 
behavior is reached only very close to the phase transition. 

The phenomenon of the suppression of the width of a non-trivial critical region is common 
in condensed matter physics 0. BCS superconductors provide the most well-known example 
of a system where criticality near T c is described by the Landau-Ginzburg mean-field theory. 



2 The whole theory of critical scaling is based on this observation. We shall also use this fact more 
explicitly when we consider a lattice theory. 
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The width, AT, of the region around T c where the mean-field description breaks down is tiny. 
There are systems where the width of the non-trivial critical region is small but measurable. 
In such systems one can observe the crossover between a mean-field and a non-trivial scaling. 

The quantitative relation between the size of the non-trivial, non-mean-field critical re- 
gion, the Ginzburg region, and certain parameters of a given system is known as the Ginzburg 
criterion. In superconductors such a parameter is a small ratio T/Ep, i.e., the width of the 
Ginzburg region is suppressed by a power of this parameter. In this paper, we show that in 
a field theory with large number of fermions, such as ([I]), the width of the Ginzburg region 
is suppressed by a power of 1/JV. 

The Ginzburg criterion can be obtained by estimating the effects of fluctuations within the 
mean-field approximation. When the fluctuations become large the mean-field approximation 
breaks down because of self-inconsistency. A measure of the importance of the fluctuations, 
or the size the corrections to the mean-field, is the value of the effective selfcoupling of the 
scalar field, A. If the coupling A is small the effect of the fluctuations is also small and the 
theory can be described by the mean-field approximation very well. However, for the 4 
theories in less than four dimensions^ fluctuations always become important close enough to 
the phase transition. This happens roughly when the coupling A^-i on the scale of m is not 
small anymore. Since \d-i has nonzero dimension equal to 5 — d, it should be compared to 
m 5 ~ d . In this way, and with the help of (|2]), one arrives at the following Ginzburg criterion 
for the applicability of the mean-field scaling: 

T 1 

N x ' 5-d V ; 

In the special case of d — 4, using (|3]) one finds: 

m>> iVlnCA/T) f ° r d = 4 - (5) 



Figure 1: An example of a graph whose contribution breaks the mean field approximation 
and the large- N expansion near T c . 

Alternatively, one can compare the size of the one-loop correction in the effective 4 
theory in d — 1 dimensions to the bare A^-i- The loop correction becomes important when 
(||) is violated.^ Indeed, consider the contribution, AA of a graph such as in Fig.|TJ to the 
effective quartic coupling A. This contribution diverges when m — > 0: 

T[A(r)] 2 

AA ~ -^FS-' (6) 

3 Here again the case of four dimensions is special. However, we are now discussing theories in d — 1 
dimensions and 1 < d — 1 < 3. 

4 In fact, the criterion (0) is well-known in the form XT c /m C 1 § as the criterion for the applicability 
of perturbation theory near a finite temperature phase transition. 
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where we integrated over fluctuations in the window between T and m with m < T. Thus, 
this one-loop contribution, AA, is negligible compared to A(T) if TA(T) <C m 5 ~ d , which 
together with (g) leads to (f|). Note also that, since A(T) ~ the contribution of 

the graph in Fig.|l] should be subleading in the large- N expansion. Therefore, the large- N 
expansion breaks down when (^) is violated. 

The Ginzburg criterion tells us that for masses m inside the window T ^> m ^> T/N x the 
mean-field scaling holds, while for smaller masses m <C T/N x (i.e., closer to the transition) 
the non-trivial d — 1 Ising scaling sets in. We see that the size of this latter, non-trivial 
critical region is suppressed at large iV.Q 

3 Lattice theory 

In this section we analyze how the effect of the suppression of the non-trivial critical region 
manifests itself in the lattice formulation of the theory. For simplicity, we shall discuss the 
case d = 3, The generalization to arbitrary d, 2 < d < 4, can be done as in the previous 
section. We consider the following lattice discretization of the theory (|]) in d = 3: 




where \ % an d X 1 are Grassmann-valued staggered fermion fields defined on the lattice sites; 
the scalar field </> is defined on the dual lattice sites, and the symbol (x, x) denotes the set 
of 8 (i.e., 2 d ) dual lattice sites x, surrounding the direct lattice site x. The fermion kinetic 
(hopping) matrix M is given by 

M x,y = \Y, V^x) [<W/i - <Wa1 > ( 8 ) 

where ?7 M (x) are the Kawamoto-Smit phases (— +^-1, The cubic lattice has L s lattice 
spacings a in spatial directions and L t lattice spacings in the temporal direction. The cutoff 
scale can be defined as A = 1/a, the temperature is given by T = 1 / (L t a) and the mass is 
m = l/(£a), where £ is the correlation length of the scalar field. To reach a continuum limit 
one has to satisfy the following two conditions: A ^> T and A ^> m. The condition A>T 
requires a lattice with sufficiently large L t 3> 1. The parameters of the action should then 
be tuned towards their critical values where the correlation length £ 3> 1. This satisfies the 
condition A ^> m. 

The condition £ 3> 1, or ma = l/£ — > 0, specifies a 2-dimensional critical surface in 
the space of 3 parameters k, A[ at and f3. One expects that a continuum limit taken at any 
generic point of this surface defines the same theory |5], |10| . (This is the meaning of the 
equivalence between Yukawa and four- fermion theories). Therefore we can fix two out of 
three bare parameters, k = Ai at = and tune a single parameter (3 to criticality: l/£ — > 0. 

5 Note that there is no proportionality constant in the Ginzburg criterion (Q). This constant would depend 
on the definition of the boundary of the mean-field region which is naturally ambiguous. The Ginzburg 
criterion tells us how this boundary moves as a function of N . 



6 



The phase diagram as a function of (3 and Ta = 1/L t looks like in Fig. |2|a. A natural 
measure of the distance from the criticality is ma = l/£. Trading the lattice parameter (5 
for ma we obtain the phase diagram of Fig. |2|b. 

classical stat. m=T 
scalar theory /' 
in d-1 

/' quantum Yukawa 
md dimensions 

ma 

Figure 2: A schematic phase diagram (a) of a lattice Yukawa theory in the plane of Ta = 1/Lt 
and the lattice action parameter f3. The solid line is the phase boundary. The dashed lines 
show the location of points where the correlation length, £, reaches L t . The same phase 
diagram (b) but the distance from the critical line is expressed in terms of a natural variable 
ma = l/£ and only one side (either symmetric or broken) is shown. Various continuum 
limits correspond to approaching the origin in (b). The slope determines the ratio T/m in 
the resulting continuum theory. The values of L t and thus of Ta are discrete, but this is of 
no importance to our discussion. 




Pc P 



The line T = m separates the regions of quantum and classical behavior or, equivalently, 
the regions of <i-dimensional and (d— l)-dimensional behavior. Below this line, when T<m, 
the correlation length £ is smaller than the extent in the time direction L t and the system 
behaves as a quantum Yukawa model in 2+1 dimensions. Above the line, when T ^> m, the 
correlation length £ is much larger than L t , the system looks like a "pancake" and behaves 
as a 2-dimensional classical statistical theory of a scalar field. 

Now consider changing the parameter (3 on a given lattice, i.e., at fixed Ta = 1/Lt, 
so that we move along a trajectory such as ABC on Fig. |3|. The effective (long distance) 
coupling A(£) follows the evolution governed by the RG equations of the Yukawa model 
from A to B. Near the point B it reaches some value, X(L t ), which, if L t is large enough, 
is given by the infrared fixed point and is 0(1/N). As we continue to increase £ from B to 
C the coupling A(£) evolves according to the RG equations of the 2-dimensional 4 theory. 
Since A (Li) ~ 1/N, in order to reach the non-trivial critical region one needs to go to the 
correlation length £ ~ L t N x , with x = 1/2 in d = 3, according to the Ginzburg criterion @. 

Thus the line of a crossover with a slope Ta/ma = N x divides the region, T > m, of a 
2-dimensional, or classical, behavior into two subregions, or windows. One window, where 
the mean-field approximation works (L t <C £ L t N x ), and another window, where the 
non-trivial critical behavior sets in (£ ^> L t N x ). 
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Figure 3: The same phase diagram as in Fig. |2|b for a Yukawa theory with large N. The 
trajectory ABC corresponds to changing some lattice parameter to approach criticality on 
a lattice with a given L t = 1/Ta. The point A corresponds to the correlation length £ = 
l/ma ~ 1, f (£) = Lt, £(D) = L t N x , and £(C) = oo. In d = 3: x = 1/2. 



4 Monte Carlo 

We performed Monte Carlo simulations of the Gross-Neveu model in d — 2 + 1 dimensions 
at finite temperature to test the results of the analysis of the previous section. We chose 
the Z2 four-Fermi model because it is relatively easy to simulate. We used a Hybrid Monte 



Carlo method described in |TT[ , which proved to be very efficient for our purposes. Since the 
chiral symmetry is discrete we were able to simulate the model directly in the chiral limit 
m = 0. This allowed particularly accurate determination of the critical properties. The 
action is that of the Yukawa lattice theory (H) with k = A = 0, and we tuned (3 to reach 
criticality. The long-wavelength (continuum limit) behavior of such a theory is determined 
by the infrared fixed point, which is the same || [K| in the more general Yukawa model (0) 
and in the Gross-Neveu model. 

We used the following two methods to optimize the performance of the Hybrid Monte 
Carlo procedure. The first method consisted of tuning the effective number of fermion flavors 
N , which is used during the integration of the equations of motion along a microcanonical 
trajectory, so as to maximize the acceptance rate of the Monte Carlo procedure for a fixed 
microcanonical time-step dr. As the lattice size was increased, the time step dr had to be 
taken smaller and the optimal N' approached N. For example, for an N = 4 theory on a 
6 x 36 2 lattice the choices dr = 0.15 and N' = 4.036 gave acceptance rates greater than 95% 
for all couplings of interest. To maintain this acceptance rate on a 6 x 60 2 lattice we used 
dr = 0.11 and N' = 4.016, while on a 6 x 80 2 lattice we used dr = 0.10 and N' = 4.012. 
In the second method the Monte Carlo procedure was optimized by choosing the trajectory 
length t at random from a Poisson distribution with the mean equal to f. This method 
of optimization, which guarantees ergodicity, was found to decrease autocorrelation times 
dramatically [13]. For most of our runs we used the average trajectory length f ~ 2.0. As 
usual, the errors were calculated by the jackknife blocking, which accounts for correlations 
in a raw data set. 

As will be seen below, we used values of the lattice coupling (3 sufficiently close to the 



8 



critical value, /3 C , so that we are close to the continuum limit A ^> m, where m is the 
thermal mass, and the scaling behavior is not affected by lattice artifacts. In addition, we 
verified that also another important physical parameter, the zero-temperature mass, m , is 
sufficiently smaller than the cutoff A. We ran on lattices with large L t , such as 20 3 with 
JV = 4 for (3 = 0.600 - 0.750, and 16 3 with JV = 24 for (3 = 0.725 - 0.875, and determined 
the value of the scaling exponent f3 m . For A^ = 12 we can use the results from [11]. 
found values in agreement with the analytical prediction (3 m = 1 + 0(1/N 2 ) for the T 



We 
= 



scaling [11 



This confirms that for our values of the coupling (3 the lattice theory remains 
in the scaling window for all range of temperatures down to T = and effects of the lattice 
are negligible. 



4.1 Exponents from finite size scaling 

The finite size scaling (FSS) analysis is a well-established tool for studying critical properties 
of phase transitions |T3|] . The critical, singular behavior in a statistical system is caused by 
the divergence of the correlation length £. On a finite lattice the correlation length is limited 
by the size of the system and, consequently, no true criticality can be observed. However, if 
the size, L s , of the lattice is large, a qualitative change in the behavior of the system occurs 
when £ ~ L s . For 1 C ^ < L s the behavior of the system is almost the same as in the bulk 
(L s = oo). However, when £ ~ L s the behavior of the system reflects the size and the shape 
of the box to which it is confined. The dependence of a given thermodynamic observable, 
P, on the size of the box, L s , is singular and, according to the FSS hypothesis, is given by: 

P(t,L.) = U*l v Q P (tL)! v ) i (9) 

where t is the distance from the critical point: t = (/3 C — (3)/(3 c ; v is the standard exponent 
of the correlation length: £ ~ t~~ v \ and Qp is a scaling function, which is not singular at zero 
argument. The exponent pp is the standard critical exponent for the quantity P: P ~ t~ pp . 
Studying the dependence on the size of the box, L s , and using @ one can determine such 
exponents. 

We simulated the model with AT = 12 fermion flavors at (3 close to the critical coupling 
(3 C . The lattice sizes ranged from L s = 12 to 40 for L t = 4, and L s = 18 to 50 for L t = 6. 
Periodic boundary conditions in the spatial directions were used. Details of the L t = 6 
runs are listed in Table [I]. To perform our study most effectively we used the histogram 
reweighting method |L4j] which enables us to calculate the observables in a region of (3 around 
the simulation coupling /? sim . 



4.1.1 Exponent v 

We used the following thermodynamic observables ]15j to determine the values of v and the 
critical coupling (3 C : 

V^4[0 3 ]-3[0 4 ], V 2 = 2[4> 2 ] - [0 4 ] , V 3 = 3[4> 2 ] - 2[0 3 ] , (10) 

where 
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One can easily find that 

V j ^(l/u)hiL a + V j (jtL 1 / v ), (12) 

for j = 1, 2, 3. At p c , i.e., t = 0, the last term in the r.h.s, Vj(0), is a constant independent 
of L s . Scanning over a range of (3's and looking for the value of (3 at which the slope of In- 
versus lnL s is j-independent, as it is in Fig. f|, we found: v = 1.00(3) and (3 C = 0.7762(15) 
for L t = 6, and v = 1.00(2) and /3 e = 0.682(2) for L t = 4. 

These values of v are in a very good agreement with the two-dimensional Ising value 
v — 1. This confirms that the behavior of the system sufficiently close to criticality is non- 
trivial (in the mean field theory: v = 1/2) and is given by the arguments of dimensional 
reduction and universality. 

4.4 
4.2 

4 
3.8 
3.6 
3.4 
3.2 

3 

2.8 3 3.2 3.4 3.6 3.8 4 

Ln(Ls) 

Figure 4: Finite size dependence of V 3 ; at (3 = 0.7762 for L t = 6. All three lines have almost 
equal slopes. 
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4.1.2 Exponents (3 m and 7 

In this subsection we consider the following two exponents: The exponent (3 m of the order 
parameter: S = {(f)) ~ t@ m , and the exponent 7 of the susceptibility: x ~ ^ 7 - According to 
eq. (^), at the critical point, t — 0, the order parameter, S, and its susceptibility, scales 
with L s as S ~ L~^ v and x ~ 

At each value of /3 we made a linear x 2_ ht f° r hi ^ and In \ versus In L s . The locations of 
the minima of x 2 /d.o.f. as a function of (3 provide estimates of (3 C . We estimated the error 
in f3 c and the critical exponents by looking at the values of (3 that increase x 2 /d.o.f. by one. 
An estimate of the error in (3 C is min(x 2 /d.o.f.) + 1, which also gives the error on the critical 
exponents. Fits at Lt — 6 for the order parameter and susceptibility gave f3 m /v = 0.12(6), 
(3 C = 0.7747(15) and j/p = 1.66(9), (3 C = 0.7750(15) respectively. Similar results at L t = 4 
are: (3 m /v = 0.16(7), (3 C = 0.6805(15), and j/v = 1.57(9), (3 C = 0.6817(15). 
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The critical exponents which we found are in a good agreement with the two-dimensional 
Ising exponents: f3 m /u = 0.125 and 7/1/ = 1.75 (in the mean field theory: (3 m /v = 1 and 
7/1/ = 2). The values of f3 c extracted from this analysis are also consistent with the estimate 
of (3 C = 0.7765(15) evaluated from the analysis of Vj. Figs. |] show the best linear fits of 
finite size dependence of ln£ and ln% at /3 = 0.7747 and (3 = 0.7753 respectively. 

-1 .86 1 1 1 1 1 1 1 5.8 1 1 1 1 1 1 




-2.02 1 ' ' ' ' ' 1 

2.8 3 3.2 3.4 3.6 3.8 4 

Ln(Ls) 



CL 




Ln(Ls) 



Figure 5: Left: the best linear fit for ln£ vs. lnL s in the minimum of \ 2 {fl = 0.7747). 
Right: the same for ln% vs. lnL s (/3 = 0.7753). 



4.2 Scaling in the broken phase 

Our task in this section is to check the analytical result of eq. (f|), which is equivalent to: 

£ < L t N x , x = 1/2. (13) 

In other words, we want to determine the position of the crossover from the mean-field (MF) 
critical behavior to the non-trivial 2d-Ising critical behavior as a function of N. 

A straightforward way to do this is to study the dependence of the order parameter, £, 
on p. Since £ vanishes at the critical point, it can be thought of as a measure of the distance 
from the criticality. We expect that for sufficiently small E, i.e., close to (3 C , this dependence 
should be given by a power-law scaling with the exponent f3 m = 1/8 of the 2d Ising model: 
£ ~ (const — /3) 1 / 8 . This corresponds to the region CD on the diagram Fig. |3|, i.e., £ 3> L t N x . 
For larger S, further away from the criticality, the MF scaling holds: S ~ (const— fi) 1 ^ 2 (with 
some other value of const). This is the region DB on Fig. ^, i.e., L t £ L t N x . For even 
larger S we should see the scaling corresponding to the fixed point of the 3d Gross-Neveu 
model fll]| : £ ~ (const — (3) 1 . In this region 1 < ^ C L t . 



Since we are interested in the boundary of the mean-field scaling region we can find £ 
from £ using a well-known relation between them: £ = l/(2£) |16[, which holds inside the 
mean-field region. Thus we avoid direct measurements of £, which are much harder than 
measurements of £. 

We studied the dependence of £ on (3 at three different values of iV = 4, 12, 24 on lattices 
with fixed L t = 6. Ideally, in order to resolve all three critical scaling regions, or windows, 
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we need to provide: 1 -C L ( C L t N x <C L s . For finite L t , L s , and N x these regions are 
squeezed, but for the values we used one can clearly resolve the MF region with the crossover 
towards the 2d-Ising region. 

In contrast to the FSS analysis where £ ~ L s , now we need to keep ^ C L s since we are 
studying bulk critical behavior. We monitored each simulation run for vacuum tunneling 
events, which signal that L s is not big enough. Away from the critical point such events were 
so rare that good measurements of £ and its susceptibility x were possible. At couplings 
near the phase transition we increased L s to suppress tunneling. Large number N, such 
as iV = 24, suppresses fluctuations and allows an accurate study within reasonable amount 
of computer resources. On the other hand, for smaller N, such as N = 4, the simulations 
were also very efficient because the crossover to the 2d-Ising behavior starts at a smaller 
correlation length. The data from these simulations is shown in Tables 0,0,f|. 

Since (3 m = 1/2 in the MF region, we fitted £ 2 with a linear function of (3. The linearity 
made the fitting procedure very efficient. We evaluated the goodness of fit for data in various 
ranges of (3 in order to find the boundaries of the MF region (see Tables ^,§,0). The drop 
in the goodness of fit when new data are added implies that the new points deviate from 
the MF behavior and they belong either to the Gross-Neveu/MF crossover region or to the 
MF/2d-Ising crossover region. Fig. ^| shows the data for £ 2 versus (3, where the straight 
lines represent the best linear fits in the MF region. All three graphs show the MF/2d-Ising 
crossovers. The graphs for iV = 12, 24 also clearly show the Gross-Neveu/MF crossover. 
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Figure 6: Order parameter squared vs. (3 for lattice theories with N = 4, 12, 24. The straight 
lines are the fits to the data in the mean-field regions. 

We used the histogram reweighting method in order to extract the value of S at the 
"border" of the MF region with the 2d-Ising region. The quality of the fit drops very 
sharply. This allowed us to use a very conservative estimate of the error, such as the width 
of the region of S in which the quality of the fit drops down to 1% (from, e.g., 90% in the 
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Figure 7: The best power-law fit (goodness 0.26) for the dependence of the boundary of the 
mean-field region, £mf> on the number of fermions, N. 

case of N — 4). The center of this region of the drop is taken as the boundary of the MF 
scaling window. We find S MF = 0.377(11), for AT = 4, S MF = 0.213(4) for N = 12, and 
Emf = 0.168(20) for N = 24. It is clear that the non-trivial 2d-Ising region is squeezed as N 
increases. In order to check the analytical prediction of eq. (§), or (|13|), we fitted our results 
to the form Emf = const • N~ x (see Fig. [j]). We found x = 0.51(3) which is in agreement 
with the analytical prediction x = 0.5. 

5 Discussion and conclusions 

In this paper we demonstrated, using analytical arguments and Monte Carlo simulations, 
that the width of the region near the finite temperature chiral phase transition where a non- 
trivial critical behavior sets in is suppressed in theories with large number of fermions, N. 
This phenomenon explains an apparent contradiction between the large- N expansion 0, |J, 
predicting mean-field scaling, and more general arguments based on dimensional reduction 
and universality [p]] , predicting the non-trivial scaling of a scalar theory in d — 1 dimensions. 

A tentative explanation of [0] was based, somewhat implicitly, on the expectation that 
the large- N expansion breaks down and cannot predict correct non-trivial exponents. Our 
key point is that we can say when, and as a function of what parameter (i.e., as a function 
of the distance from the criticality, m/T) the large- N expansion breaks down. What is, 
perhaps, even more important, is that we show that there does exist a region where the 
large- N expansion is valid. This region (in the space of m/T) of the mean- field behavior 
squeezes out the region of the true non-trivial critical behavior.^ 

One can also look at the whole problem as a question of the order of limits. If the limit 
iV — > oo is taken before m/T — > 0, the non-trivial region disappears, and all critical behavior 

6 It was conjectured in |l7j that the Ginzburg region of non-trivial scaling could turn out to be small. 
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is given by mean-field theory. If, on the other hand, the limit m/T — > is taken at fixed N, 
the true non-trivial scaling will hold. These effects can be clearly seen in our Monte Carlo 
studies. 

A helpful analogy can be drawn with the question of the long distance behavior of the 
scalar field correlator in the original 2-dimensional Gross-Neveu model with a 2/(1) symmetry 



at large N [18|. The large- N expansion predicts spontaneous symmetry breaking and long- 



range order, which is in clear contradiction with the Mermin- Wagner theorem. As was shown 



by Witten [Tj|, there is an "almost long-range order" in the system, i.e., the correlator falls 
off like r~ l / N . Looking at this expression, one can easily see the interplay between the limits 
r — > oo and iV — > oo. 

We applied our arguments to a specific example of a Yukawa theory, but the mechanism 
responsible for this phenomenon is clearly more universal and may apply, in particular, to 
QCD, provided that the number of colors N c is large. The Yukawa, or the Nambu-Jona- 
Lasinio, theory is well-known to provide a very good description for the phenomenon of chiral 
symmetry breaking and restoration. One can then think of the theory we considered as an 
effective description of the degrees of freedom participating in chiral symmetry breaking. 

As we have seen, the role of the fermions is to screen the effective self-coupling of the 
scalar field, A. The strength of this effect depends on two factors: (i) large N, and (ii) 
large window of scales between the cutoff of the effective theory, A, and the temperature, T. 
Let us see if these conditions are fulfilled in QCD. The scale of the spontaneous symmetry 
breaking is of order A ~ 1 GeV, while T c pa 160 MeV. Thus, there is almost an order of 
magnitude window between A and T, which is presumably sufficient to drive the effective 
self-coupling of the scalar field to its infrared fixed point at the scale of T c . How small this 
value is now depends on the number of the fermions (the condition (i)). In QCD the value 
N c = 3, though not very large, can be considered large in some cases. It would be interesting 
to see if this phenomenon could be rigorously shown to occur in the limit N c — ► oo, which is 
very plausible. 

The effect of the suppression of the width of the non-trivial critical region in QCD may 
lead to the following prediction: only mean-field (but no non-trivial 0(4)) scaling behavior 
could be seen because of non-zero quark masses, m q . On the one hand, the mean-field 
scaling would hold until relatively large thermal correlation length. On the other hand, this 
correlation length in the real world is limited by the quark masses, m q . 

The following crude estimates can serve to illustrate this point .[] A true criticality is 
never reached near T c , because m q plays the same role as the external ordering magnetic 
field in a ferromagnet. The largest thermal correlation length in units of 1/T, T/m, can be 
estimated using the analogy to a ferromagnet and the mean-field value of u/({35) = 1/3, asQ 

160 MeV \ 1/3 




~ 3. (14) 



5 MeV J 

The fact that this number is not large could be guessed by observing that the zero tempera- 



7 To make these estimates quantitative one needs to find numerical coefficients, which are determined by 
non-universal dynamical properties of QCD and by the definition of the boundary of the mean-field region. 

8 We use the definition of the exponents: m ~ t u , E ~ t@ , and E ~ m^ S , together with the central 
postulate of the scaling theory: critical properties are determined by the correlation length, 1/m. 
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ture pion masses, which are driven by m g , are as large as T c . This largest correlation length 
may turn out to be smaller than the one required for the crossover to non-trivial scaling 
region: T/m ~ iV c ln(A/T c ) ~ 6, according to ([5|). In this case, the (near-)critical behav- 
ior observed in the window allowed by non-zero quark masses (according to fll4|), roughly: 
1 < T/m < 3) will be given entirely by the mean-field scaling. 
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Table 1: Simulations for the FSS analysis with L t = 6 and iV = 12. 



L s Aam trajectories 

18 0.7744 130,000 

18 0.7763 100,000 

24 0.7744 100,000 

24 0.7764 80,000 



L s /9sim trajectories 

30 0.7744 120,000 

30 0.7764 130,000 

40 0.7754 140,000 

50 0.7754 152,000 



Table 2: The data for the scaling of E in the broken phase. N — 4. 



p 


E 


L s 


trajectories 


P 


E 


L s 


trajectories 


0.595 


0.5008(4) 


36 


30,000 


0.635 


0.3736(7) 


36 


70,000 


0.600 


0.4871(4) 


36 


40,000 


0.640 


0.3524(5) 


36 


110,000 


0.605 


0.4724(4) 


36 


40,000 


0.645 


0.3297(8) 


36 


130,000 


0.610 


0.4578(4) 


36 


40,000 


0.6475 


0.3162(9) 


36 


150,000 


0.615 


0.4426(4) 


36 


50,000 


0.650 


0.3016(10) 


48 


54,000 


0.620 


0.4268(4) 


36 


54,000 


0.652 


0.2858(15) 


60 


41,000 


0.625 


0.4106(3) 


36 


70,000 


0.654 


0.2647(26) 


80 


25,400 


0.630 


0.3928(4) 


36 


64,000 


0.656 


0.2398(30) 


80 


30,000 



Table 3: The data for the scaling of E in the broken phase. N = 12. 



p 


E 


L s 


P 


E 


L s 


P 


E 


L s 


0.640 


0.5112(2) 


30 


0.710 


0.3608(4) 


30 


0.764 


0.2048(2) 


60 


0.650 


0.4911(3) 


30 


0.720 


0.3368(4) 


30 


0.766 


0.1951(3) 


60 


0.660 


0.4701(2) 


30 


0.730 


0.3119(4) 


30 


0.768 


0.1844(3) 


60 


0.670 


0.4490(3) 


30 


0.740 


0.2847(4) 


30 


0.770 


0.1725(3) 


60 


0.680 


0.4276(2) 


30 


0.750 


0.2557(8) 


30 


0.772 


0.1575(5) 


60 


0.690 


0.4061(3) 


30 


0.755 


0.2393(5) 


30 


0.774 


0.1388(7) 


60 


0.700 


0.3836(2) 


30 


0.760 


0.2199(5) 


30 
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Table 4: The data for the scaling of £ in the broken phase. N = 24. 



P 




T 

L s 


trajectories 


p 


Zj 




trajectories 


0.6700 


0.4756(4) 


24 


3,000 


0.7700 


0.2617(4) 


24 


8,000 


0.6800 


0.4560(4) 


24 


3,000 


0.7800 


0.2345(5) 


24 


8,000 


r\ r* c\r\r\ 

0.6900 


0.4359(4) 


24 


3,000 


0.7900 


0.2035(5) 


24 


10,000 


0.7000 


0.4162(4) 


24 


3,000 


0.7950 


0.1874(3) 


36 


15,000 


0.7100 


0.3963(4) 


24 


3,000 


0.8000 


0.1676(3) 


36 


20,000 


0.7200 


0.3760(5) 


24 


3,000 


0.8025 


0.1568(4) 


36 


20,000 


0.7300 


0.3551(5) 


24 


3,000 


0.8050 


0.1437(4) 


48 


24,000 


0.7400 


0.3325(5) 


24 


3,000 


0.8075 


0.1267(10) 


60 


12,000 


0.7500 


0.3101(5) 


24 


3,000 


0.8100 


0.1074(15) 


60 


22,000 


0.7600 


0.2866(3) 


24 


8,000 


0.8112 


0.0914(28) 


72 


10,000 



Table 5: Goodness of linear fits of S 2 vs. (3 for various ranges of /3. N = 4. 



no. of points 


(3 range 


goodness 


no. of points 


(3 range 


goodness 


4 


0.595 - 0.610 


0.92 


8 


0.595 - 0.630 


0.96 


5 


0.595 - 0.615 


0.95 


9 


0.595 - 0.635 


0.21 


6 


0.595 - 0.620 


0.98 


10 


0.595 - 0.640 


10" 5 


7 


0.595 - 0.625 


0.99 


11 


0.595 - 0.645 


io- n 




Table 6 


: The same 


as Table § but for N = 12. 




no. of points 


(3 range 


goodness 


no. of points 


(3 range 


goodness 


5 


0.70 - 0.74 


4 x lO" 3 


6 


0.72 - 0.760 


0.60 


6 


0.70 - 0.75 


5 x 10~ 3 


7 


0.72 - 0.762 


0.40 


5 


0.71 - 0.75 


0.20 


8 


0.72 - 0.764 


5 x 10~ 4 


4 


0.72 - 0.75 


0.67 


9 


0.72 - 0.766 


2 x 10~ 15 


5 


0.72 - 0.755 


0.77 









Table 7: The same as Table [| but for iV = 24. 
no. of points (3 range goodness no. of points (3 range goodness 



7 


0.74 - 


0.795 


6 x 10" 4 


5 


0.76 - 


0.795 


0.70 


8 


0.74 - 


0.800 


1.2 x 10~ 3 


6 


0.76 - 


0.800 


0.38 


6 


0.75 - 


0.795 


0.28 


7 


0.76 - 


0.8025 


0.17 


7 


0.75 - 


0.800 


0.24 


8 


0.76 - 


0.8050 


1.2 x 10~ 8 


3 


0.76 - 


0.780 


0.68 


9 


0.76 - 


0.8075 


5 x 10~ 16 


4 


0.76 - 


0.790 


0.91 
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